Dislocation jamming and Andrade creep 



(N 
O 

o 

(N 



43 

o 

CD 



M. -Carmen Miguel, 1 Alessandro Vespignani, 2 Michael Zaiser, 3 and Stefano Zapperi 4 

1 Departament de Fisica Fonamental, Facultat de Fisica, 
Universitat de Barcelona, Av. Diagonal 647, E-08028, Barcelona, Spain 
2 The Abdus Salam International Centre for Theoretical Physics, P.O. Box 586, 34100 Trieste, Italy 
3 Center for Materials Science and Engineering, University of Edinburgh, 
King's Buildings, Sanderson Building, Edinburgh EH93JL, UK 

4 INFM UdR Roma 1 and SMC, Dipartimento di Fisica, 
Universita "La Sapienza" , P.le A. Moro 2, 00185 Roma, Italy 
(Dated: February 1, 2008) 

We simulate the glide motion of an assembly of interacting dislocations under the action of an 
external shear stress and show that the associated plastic creep relaxation follows Andrade's law. 
Our results indicate that Andrade creep in plastically deforming crystals involves the correlated 
motion of dislocation structures near a dynamic transition separating a flowing from a jammed 
phase. Simulations in presence of dislocation multiplication and noise confirm the robustness of this 
finding and highlight the importance of metastable structure formation for the relaxation process. 
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Andrade reported in 1910 that the creep deformation 
of soft metals at constant temperature and stress grows 
in time according to a power law with exponent 1/3, i.e. 
7 ~ i 1 / 3 where 7 is the global strain of the material 0. 
More generally, the creep deformation curve usually fol- 
lows the relation 7 = 70 + /Si 1 / 3 + ret, where 70 is the 
instantaneous plastic strain, /Ji 1 / 3 is Andrade creep, and 
nt is referred to as linear creep j2| This same be- 
havior has since been observed in many materials with 
rather different structures leading to the conclusion that 
this should be a process determined by quite general prin- 
ciples, independent of most material specific properties. 
In crystalline materials, Andrade's and linear creep find 
their microscopic origin in the dynamics of dislocations, 
the topological defects responsible for their plastic defor- 
mation ||, ||, Q]. Plastic flow only occurs when the ex- 
ternally applied stress overcomes a threshold value, the 
yield stress of the material, such that large-scale disloca- 
tion motion may take place. Despite various arguments 
proposed within the dislocation literature |^|, |^, [| ^, [?], 
there is no general consensus on the basic mechanism to 
explain Andrade's law. Mott || attributed the power law 
to an athermal cooperative process taking place close to 
the yield stress, but his idea was not worked out. Later 
explanations have always focused on thermally activated 
processes |6| [?J and rely on assumptions that may not be 
fully warranted M. 

Under the action of external stress, dislocations tend to 
glide cooperatively due to their mutual long-range elastic 
interactions, which can be attractive or repulsive. The 
anisotropic character of dislocation-dislocation interac- 
tions contributes at the same time to the formation of 
spatial dislocation structures observed in transmission 
electron micrographs of plastically deformed metals || . 
As a result of these peculiar interactions, self-induced 
constraints build up in the system and the motion of dis- 
locations may eventually cease. Small variations of the 
external loading, the density, or the dislocation distribu- 



tion can, however, enhance dislocation motion in a dis- 
continuous manner. This gives rise to rather complex and 
heterogeneous spatio-temporal patterns of plastic flow, 
which have been observed experimentally in the form of 
slip lines and slip bands emerging on the surface of met- 
als jll| ElL fill or in the acoustic emission activity of ice 
crystals |13[|. Other processes commonly taking place in 
plastically deforming crystals such as hardening, fatigue, 
or plastic instabilities |l4|, ^5|, are further consequences 
of the dislocations intriguing behavior. 

In this Letter, we study the temporal relaxation of a 
relatively simple dislocation dynamics model through nu- 
merical simulation. In particular, we consider the behav- 
ior of parallel straight edge dislocations moving in a sin- 
gle slip system under the action of constant stress. We 
show that the model gives rise to Andrade-like creep at 
short and intermediate times for a wide range of applied 
stresses, without invoking thermally activated processes. 
At larger times the strain-rate, which is proportional to 
the density of mobile dislocations, crosses over to a linear 
creep regime (steady rate of deformation) , whenever the 
applied stress is larger than a critical threshold er c , or de- 
cays exponentially to zero. These results suggest that a 
possible interpretation of the creep laws could be found 
within the general scenario wrapping a dynamic phase 
transition from a flowing to a jammed dislocation phase. 
The "jamming" scenario has been recently proposed |ll| 
to understand a broad class of non-equilibrium physical 
systems (granular media, colloids, supercooled liquids, 
foams) which, in spite of their differences, exhibit com- 
mon properties such as slow dynamics and scaling fea- 
tures near the jamming threshold. When jammed, these 
systems are unable to explore phase space, but they can 
be unjammed by changing the stress, the density, or the 
temperature. To further explore the analogies of disloca- 
tion motion and these so-called jammed systems, we con- 
sider the influence of dislocation multiplication, due for 
instance to the activation of Frank- Read sources [|[ ||[ |J] 
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during the deformation process, and thermaUike fluctu- 
ations on the dynamics of the dislocations. Dislocation 
multiplication favors the rearrangements of the system 
and induces a linear creep regime (flowing phase) at lower 
stress values, but it does not affect the initial power-law 
creep. The introduction of a finite effective temperature 
T has a similar effect. 

We consider a two-dimensional (2d) model represent- 
ing a cross section of a single-slip oriented crystal where 
N point-like edge dislocations glide in the xy plane 
along directions parallel to the x axis. Dislocations 
with positive and negative Burgers vectors (the topolog- 
ical charge characterizing a dislocation) b„ = ±bx are 
assumed to be present in equal numbers, and the ini- 
tial number of dislocations is the same in every realiza- 
tion. Several 2<i-models containing similar basic ingredi- 
ents have been proposed in the literature in the last few 
years [[l3| [l?], [T^, |l9) . An important feature common to 
these models is that dislocations interact with each other 
through the long-range elastic stress field they produce 
in the host material. An edge dislocation with Burgers 
vector bx located at the origin gives rise to a shear stress 
a s at a point r = [x, y) of the form |j, Q 



a s (r) = D 



x(x 2 - y 2 ) 



(x 2 + y 2 ) 2 



(1) 



where D = h/2 , k(\ — v) is a coefficient involving the shear 
modulus /j, and the Poisson ratio v of the material. We 
further assume an overdamped dynamics in which the 
dislocation velocities are linearly proportional to the local 
forces. Accordingly, the velocity of the nth dislocation 
along the glide direction, if an external shear stress a is 
also applied, is given by 



= b n (y~] 0- s (r nm ) + a), 



(2) 



where Xd is the effective mobility of the dislocations pQ] 
(Xd 1 /b the effective friction per unit length) and r nm = 
r n — r m the relative position vector of dislocations n and 
to. Periodic boundary conditions are imposed in the di- 
rection of motion (i.e. the x axis). In order to take 
correctly into account the long range nature of the elas- 
tic interactions (Eq. [IJ) , we sum the stress over an infinite 
number of images. This sum can be performed exactly 
and the results are reported in Ref. jjj. When the dis- 
tance between two dislocations is of the order of a few 
Burgers vectors, linear elasticity theory (i.e. Eq .(]]])) 
breaks down. In these instances, phenomenological non- 
linear reactions, such as the annihilation of a pair of dis- 
locations, describe more accurately the real behavior of 
dislocations in a crystal [jlj| [l?], |l8[ fL9| . In our model, 
we annihilate a pair of dislocations with opposite Burgers 
vectors when the distance between them is shorter than 
y e (for Cu, y c ps 1 .6nm plj). In the following, we measure 
all lengths in units of y e , time in units of t — y 2 / (\Db 3 ), 
and stresses in units of Db/y c . 
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FIG. 1: The strain rate relaxation for different applied stresses 
at T = for a system of size L — 100y c . The initial density 
of dislocations is around 1%. The solid line is the best linear 
fit of the a = 0.01 curve and yields drf/dt ~ t~ ' 69 . 



To analyze creep relaxation, we integrate numerically 
the N coupled equations using an adaptive step size fifth- 
order Runge-Kutta algorithm. Simulations start from a 
configuration of Nq dislocations randomly placed on a 
square cell of size L. We have considered two different 
box sizes L = 100z/ c , and L = 300y c , with initial numbers 
of dislocations No — 400, and iVo = 1500, respectively. 
We first relax the system until it reaches a metastable 
arrangement, where the remaining dislocation density is 
around 1%. Next we apply an external shear stress and 
let the system evolve. The results are typically averaged 
over several (100-400) random initial configurations. 

In Fig. |l| we report the plastic strain rate of the ma- 
terial, defined as dj/dt ~ ^j^biVi, for different values of 
the applied stress. The strain rate decays as a power law, 
with an exponent close to 2/3, in agreement with An- 
drade's law. For high stresses the power law relaxation is 
followed by a plateau indicating the onset of a linear creep 
regime. The crossover time increases as the stress de- 
creases, and for stresses lower than a ~ 0.0075 — 0.01 the 
plateau disappears and the creep decays exponentially 
to zero. In Fig. || we display the steady-state strain rate 
as a function of stress. This plot suggests the presence 
of a non-equilibrium phase transition between a moving 
and a jammed stationary state controlled by the applied 
stress. Using a c = 0.0075, we find dj/dt ~ (a — a c )P 
with /3= 1.8 ±0.1. 

From our model and the analysis of our data, we pro- 
pose a new explanation of Andrade creep. In the course 
of time, most dislocations tend to arrange themselves into 
metastable structures. These structures often consist 
of small-angle dislocation boundaries separating slightly 
misoriented crystalline blocks, whose stress fields are 
screened on large length-scales, of dislocation dipoles, or 
of far more complex dislocation arrangements fOjfl . The 
stress field generated by a fraction of the dislocations in 
these metastable configurations conserves, however, the 
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FIG. 2: The steady- state value of the strain rate as a function 
of the applied stress. The inset shows that the strain rate 
scales as (<r — Uc) 13 , with (3 = 1.8. 



initial long-range character, and this enforces the contin- 
uation of the relaxation process. If the applied stress is 
close to the critical threshold, dislocations rearrange in 
a correlated and intermittent manner p3[ exploring fur- 
ther configurations, yielding Andrade law as the average 
result. Above the critical stress, the system eventually 
enters a linear creep regime in which dislocation struc- 
tures are coherently moving with a velocity that grows 
with the applied stress. On the contrary, for stress val- 
ues below the critical threshold, dislocations cannot get 
around the self-induced constraints built up by their mu- 
tual (attractive/repulsive) interactions and thus they are 
unable to further explore configuration space, the signa- 
ture of a jammed system The closer is the applied 
stress to the threshold, the longer is the extent of the 
power-law regime before the system falls either in the 
flowing or in the jammed state. Precisely at the critical 
point and in an infinite system, the Andrade power-law 
would last indefinitely. 

Within the proposed scenario, it appears natural to 
ask whether two important effects, namely dislocation 
multiplication and thermal fluctuations, affect this be- 
havior in a relevant way. During plastic deformation, 
in fact, new dislocations are created within the crys- 
tal. It is widely believed that the Frank-Read multiplica- 
tion mechanism || |[ Q is the most relevant for disloca- 
tions gliding under creep deformation. Since Frank-Read 
sources cannot be directly simulated in a two dimensional 
model, we employ a phenomenological procedure, intro- 
ducing dislocations pairs with a rate r proportional to the 
external stress. Similar multiplication mechanisms have 
been successfully used in the past []l3[ |l7|, |l8|, ^9|. Ther- 
mal fluctuations are accounted for by adding a Gaussian 
random force per unit length rj n (t)/b to Eq. |2|. The force 
has zero mean and its correlations are given by 



{Vn(t)Vm(t')) = k B Tx^S n , m S(t - t'), 



(3) 



where T is an effective temperature characterizing the 
strength of the fluctuations 1 22 . This random force could 



FIG. 3: The effect of the dislocation multiplication rate r and 
of the effective temperature T on the plastic strain in a system 
of size L = 100y e with an applied stress a — 0.0075. Andrade 
creep (7 ~ t 1 ^ 3 ) crosses over to linear creep (7 ~ t) at large 
times, with a crossover depending on r and T. 



also mimic, as a first approximation, the influence of dis- 
location motion in other slip systems that may be active 
simultaneously in the material (see j^, p3[). 

We have performed a series of simulations for differ- 
ent multiplication rates and effective temperatures (mea- 
sured in units of Db 3 /fee) when the applied stress is close 
to the critical value (i.e. a = 0.0075). We used the 
same procedure and initial conditions described previ- 
ously and, given the additional fluctuations present in 
the strain rate, we focus our attention on the integrated 
plastic strain 7. As shown in Fig. Q a linear creep regime 
appears after a crossover time which decreases with r 
and T. Nevertheless Andrade creep (i.e. 7 ~ i 1 / 3 ) still 
persists at shorter times. A visual inspection of the dis- 
location arrangements during the deformation is useful 
to understand the origin of the various creep regimes. 
At low T, we observe the presence of slowly relaxing 
metastable structures as for T = 0. On the other hand, 
at much higher effective temperature (i.e. T ~ 1) these 
structures have completely broken up (melt) , the flow of 
dislocations becomes fluid-like, and Andrade creep disap- 
pears. In order to substantiate quantitatively this state- 
ment, we measure the angular correlation between dis- 
locations. For each dislocation pair with coordinates fj 
and fj , we define 9 as the azimuthal angle with respect to 
the y axis of the vector fij (i.e. 9 = arccos(j/ • ?%/|t%|))- 
Thus two dislocations placed in the same wall have ei- 
ther 9 ~ or 9 ~ 7T, a dislocation dipole is characterized 
by 9 ~ tt/4, 37r/4 and two dislocations in a pileup yield 
6 ~ tt/2 @, |, f§. The distribution P(9) is obtained after 
averaging over all dislocation pairs in several realizations 
(~ 100) of the dynamics at a given instant. Fig. [| shows 
P(9) for different values of T. At low effective temper- 
atures we observe five roughly equivalent peaks, corre- 
sponding to walls, dipoles and pileup configurations. As 
the effective temperature increases, the peaks at 9 = 0, tt 
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FIG. 4: The angle distribution as a function of the tempera- 
ture at a = 0.0075. 



and 9 ~ 7r/4,37r/4, corresponding to walls and dipoles, 
progressively disappear, while the 6 = n/2 peak remains. 
This peak indicates that even in the high temperature 
regime the dislocations are correlated in the x direction, 
the only possible direction of motion. Apart from this no 
other structure remains in the system. 

In conclusion, our results capture several features of 
the phenomenology observed in creep deformation exper- 
iments. Andrade's scaling and the creep curve appear 
to be controlled by a non-equilibrium phase transition 
between a jammed and a flowing dislocation state. De- 
termining the critical stress value, and characterizing its 
behavior (possible dependence on density and temper- 
ature) is of utmost importance for practical purposes, 
since it establishes the mechanical strength of crystalline 
materials. 
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